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ABSTRACT 
, The use of parameter plane curves of output voltage magnitude versus 
one of the parameters as a means of obtaining the frequency response of 
nonlinear networks was investigated for networks containing one nonlinear y 


element, which was a function of the magnitude of the voltage across it. 


~ 


The results obtained were checked by digital simulation and it may be 
concluded that the method is reliable. The method was also presented 
as a way of obtaining information for the reshaping of network frequency 


response by introduction of the proper nonlinear elements, 


| 
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ie Introduction to Parameter Plane Techniques 

Parameter plane work was first developed to investigate the effects 
of the variation of two elements (parameters) upon the location of the 
roots of the characteristic equation of control systems. 

Because of the nature of the characteristic equation, this work was 
later extended to the analysis and design of networks with linear elements, 
especially because of the possibility of controlling and modifying the 
characteristics of the network by the variation of two of its elements. 

The referred parameters are considered to be varied slowly and to 
be time invariant, as far as the development of parameter plane techniques 
are concerned; but the present work extends one of the aspects of this 
utilization, by using them in the study of the frequency response of net- 
works with nonlinear elements, with time variations of their own. 

Basically the parameter plane method can be applied to the analysis 
of any polynomial with variable coefficients. 


Consider the equation 
fos) = on hey Ss te (1-1) 
Keo 


where ., are real coefficients and s is the complex frequency variable 
(Ss]=2 {wy 


Considering the Chebyshev functions of the first and second kind, 


RCH=C) i () = 


and 


UceH=e) U.@) a) 


and because g = - Wy + { U Vi- 2% , it follows that: 


rw) 


ee a (hes) +y ie Uc) 9 





i. and tse obey to the recursion relationships given by the tables of 


the Chebyshev functions 


Table I-1 


Ter ath, @) 7,2) 


Uenteh 2 TU) = UG) 
i, (2) =! i (t)=% 


U, (2) = 0 U @&) = / 


Equation 1-4 results from the fact that: 


sa - mh & +f 4 /, - 2 


where WwW, is the undamped natural frequency and < is the damping co- 


efficient; je 
also Sz w, (om + {sim #) = Oy, L > “4 $= és @4) 


k , JK? K _— 
Then S = Wy & = Wy, (on Ke + } sim Ko) ; 


an (Kk en”'& E)) 


defining: T, G2) = cr K+ 


Ulcer) = sie ke 2 stu (Kon (C8) 
Sim Sim (cas “TE z)) 





the result is: 


; se we LC®) + Vi-8* v.62) (1-4) 


Using the recursion relationships for the i and CU functions: 


sous (ey 7 ie ei)” U, (2)] nes 


Using these expressions for s in one original equation for the poly- 
nomial in s, there results a complex expression, for which the real and 
imaginary parts are equated to zero independently, resulting in two 


equations 


et 


z. a, un (-1) T, (5) = 0 


(1-6) 
a va K+ 
2. A. Wy (1) gb (z)=0 | (1-7) 


Using Table I-1l, one obtains the relationship 


T(t) = t UG) - U., @) (1-8) 
2 and finally for (1-6) and (1-7) 
s — K ; 
f Z ti) a, Ww" i i (4) = 0 (1-9) 
K= 0 
— K 
Z. ape Uepse —~ 
K=0 


2 





These expressions lead to the final expression of ™ and (3 in terms 


of Z and ~,, for the cases when the coefficients “, are of the form 


Bie = b,. tee. 8 + Ay y » where s- and | 


show up only in linear combination. 

For a reasonably general situation, the coefficients “x are 
considered to be functions of the two variable parameters, with the 
form 

b % i (? ie Be fp + fx 

when 6, Cy, Ax» fe are constant and real doettieionts and ~° and 
fp the two real variable parameters. ; 

Using the previous derivations for ie and the te and Ui. tunetions, 


uw 


the following expression is obtained: 


fis) = x a. wn T.¢t)+ jVi-F pa Ay Wr Ud aay 


Manipulation of the relations in Table I-1' leads to equations 1-6 
and 1-7, and by substitution of the general expression for , into them, 


there results, 


i 
G 


Bw +, (> + D, x 3 ws f, 


(4-19) 


W 
© 


Bx +0, 0+ eps Fy 


(1-13) 


These were obtained from the equations 


x ©) b w"| Tz) + ed es) eye TS) « (1-14) 
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and 
Ww 


«EC go Y)+ AZ Eye wo VG) + 
— ro 


(1-15) 


a x (3 iC AG Me i UC) +2 fw" U,@) =O 
Keo = 


and by calling 


Bj 22 ul « b, w* 7%) 
Bz - me i) b, ws UG) 
Keo 


C= L en", ot EE) 


: C22 bie, ut U, @) 
K=O 
” K 

: D = (1) dy ws Lo (1-16) 
kro 


By elimination of © , an expression for (3 is obtained in the fol- 


lowing fashion, from (1-12) and (1-13): 


shat + Gc G a Fo t C2 2 or 
B,+D, B2 + Dz 7 


1l 





Aunt « (Agr Ag)¢ Aa <2 


where 


Ax. 2 CD ORs 9 2,28 Bite. -- 


Also 


he = ~CApp+ Oc) ¢ IC Arp + Acs) -4 Ben Ary 


2 Be» 


Ka 


and by defining 


one obtains: 





(C? sa. [- (Arp + Arg) # Re | | 


In a similar way, eliminate - and obtain 


where 


Kee = (Ar. i Cee 4 An Lg 


It can be shown (Hollister, Ref. 1) that a = i ‘ 


the common value and obtain 


) 


‘ra ; (Aro A..)* RI 





[ = Ogg? fg) 4 A., Ar 


etc., 


Call / 


(1-17) 


(1-18) 


(1-19) 


(1-20) 


Cw 


(1-21) 





Using the fact that Bias = - os ,» and from equations (1-20) 
and (1-21), which represent four possible pairs of solutions, Hollister 
c shows that the only ones which are solutions for the original (1-12) 


‘and.-(1-13) equations in ~< and leg » are 





2 Bap L 


Pan - (1-214) 
Bx’ ’ oie [= (eo A...) qi!) e| 
where kK = 0 andk = 1 sive the two possible simultaneous solutions for 
eq. (1-12) ‘and (1-13). : 
These ape the, two basic-equations that were used to develop digital 
* computer programs on use in the parameter plane. By fixing ~% or wy» 
one can foive the pquneten ERE e and fs for particular values of 
“Wn OF 3 respectively, and dbtain families of curves of & versus 
= - » used in parameter plane analyses. This basically consists in 
‘seeing. ftom these curves which relations should be maintained between 
e and e to locate the roots of’ the polynomial in places correspond- 


ing to desired undamped frequencies ¢ Wr = We ) and desired damping 


. 


ratios (-& = % y. 





Dea: Curves with Magnitude as One of the Coordinates 


In the present work primary interest is in the curves relating the 
magnitude of a given transfer function and one of the parameters in the 
system, for successive values of real frequency W . 

This is obtained as follows (Ref. 1); consider the transfer function 


for a given physical system 


<=} 


G(s) = Gn Spt Gaal + se SH He (2-1) 


=] 
An $"+ i p tre t aS 4 Ao 





where coefficients are of the form 


: a = Dp eo Ed ao +h (2-2) 


dr = usb Bre npr ty (2-3) 


Then 


ee 


Gis) =e (te tte B+ Be 2B + Ae) 5” 
da (a,% i. bp G+ Cre A + dp) s? 


To study the real frequency variations, put S = in the transfer func- 





fa 


tion and obtain 


wen 
MA 


: 
Gijnye dee ON [wanes fe Ot fae + he) tres, 
= wate A ae bp 3 + crx > oo Ap) + 


po 


yale!) 


ve® ) “ 1 X+ fy Bt Fu%A+ &., ) (2-4) 
‘ Vo ( aii 
pl > ey” a w (ap% + bp Arce +dp) ¢ 
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For each particular value of () chosen, @)= CU, ), the following 


definitions are made: 
cen 


An aD it ‘Apap De ape | 


and sata for CA ; Dr Q 


rae 
A; Sut WU, ap 5 Be 


and hee for An 5 D: . 


ne 


ie 
s 
s 


even 


Nhe Ate R : 

i 4 > (-1)* Wy ER 2 and SIMiLAnuyY FOR The Gy ! Hn 
R20 

and 
odd 


(Me d(k-!) gen 
E. - (-!) Wy, Ep and Sranied Rey FOR ' eo, e 


Finally there results: 


(jos) _ jek +BFi+xp Tat | + gl: “phi +x +H] (e-5) 


I Tanta In FXG, + Dy] 4g [Alt Bi BCE FDL T 


Per forming G(4He) * G(-4.) one obtains the magnitude of the 


transfer function defined by i 


M= Go) 


es 


2 2, 
Me (xEnt ah +B Gath) bE eB + XRG: +H) 
© (KAR + BB +48 Cnt On)* + (AL HBB FA BC; HDI? (2-6) 


These manipulations provide an equation in & and A for M?, which 


has the general form 


@ (4:8) -[M"-6(«,~)] = ° 


and by expanding each of the components of (2-6) and isolating 6 (as 











a chosen variable parameter) there results: 


(M*R, ~Fa)arell “W, - W2)2 +MV-V, = 0 (2-7) 
where R= (Bn + Bi) + Ix (Br Cnt Bele) 4 aC + Ci ‘ 
Ree (BHF) + aKlhon FG) rA(G + ) 


W, = 2(Bn Dy + Di Be)+2alby By tr Be +C2 Dé #Cx Dn) + 
+2 a* (AnCn t Accs] 


Wa. 20 Hye ibd) + 2 «(EA HELE HOM HM 
7 pa*(En Cy, +bG,)| 


V, = CODE) 4 a(n Oye) aif 


V2 . (pot H;'] t Joa (Enty rE AJP E8| 








As the final step, solve the quadratic equation for Ga and obtain 


a= (4 wi-w) + VMtwi- Wa) 4 (28, -Re)(miy 4) (2-8) 
4,2 VA (M7R, ~R ) 

As was stated earlier, equations 1-20A and 1-21A could be used to 
write digital computer programs for parameter plane analysis. The same 
may be said of equation 2-8, which provides the possibility of obtaining 
a family of constant frequency curves of voltage magnitude versus one 
of the original parameters “ or 3 » Or any combination of the four 
variables ( « » »Mandw ) by fixing one of them and combining the 
other three in any desired way. 

This was in fact the work of Rommel and Cradit (Reference 4) which 
prepared computer programs designed to give the possible different plots 
resulting from the combination of the four variables referred to above, 
and Glavis (Reference 3) who perfected these programs and combined them 
into one unique deck prepared to give any combination of parameters or 
variables which may be of interest (See Appendix I). 

This program, which was used in the preparation of this thesis, is 
included as Appendix I PARAMS (Digital Computer Program). Using it, 
frequency curves are obtained for magnitude as function of one of the 
parameters, maintaining the other parameter fixed; this corresponds to 


the PARAM 5 Part of the program cited above. 


7. 
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“3.1. Description of the Method. . 


3, Application of Parameter Plane Curves to Frequency Response. 


Having already available a computer program to solve the system 
cnanpter tineciot ant determine ike iagat thde~ebesue-frequency response 
did not solve all the questions about obtaining plots Ser analysis of° 
systems from the respective transfer functions, because all the programs 
were prepared for work on cases where all the elements are time invariant. 

The task of —, the response of a system with nonlinear ele- 
ments, or even etihinvariant elements, was not easy to’ solve analytically. 
Therefore, it was thought that the parameter plane curves of the PARAM 
5 type (voltage magnitude versus one of the parameters, with ™“ varying 
from curve to curve) could be used to obtain quickly the frequency response, 
by graphically reading the valitee of voltage magnitude corresponding to 
the particular value of the parameter where it intersects each frequency 
curve, ‘and plotting a new graph of magnitude versus frequency, for that 
particular value of the parameter. . 

This procedure also makes it very practical to use the graphical 
representation of the nonlinearity if it happens to be a function of the 
magnitude, and by use of superposition of graphs, to obtain the frequency 
response for the nonlinear cases, 

The purpose of this werk is kd show how this is done, by means of 
some examples for the case of electrical networks, and to prove that the 
aetked is good in the solution of this ‘type of problem, giving results 
which may be considered within a accuracy. 


. 


' The method described here to obtain the ‘frequency response of electri- 


cal networks containing one nonlinear element, may be summarized as fol- 


lows: 


1. Obtain the transfer function of the network, in terms of 
output function ie a function of the input. The most direct applica- 
tion of the method occurs when the nonlinearity present is a function of 
the magnitude of the voltage across it (for example, temperature-sensi- 
tive resistors or voltage-sensitive capacitors). This voltage is the 
output function wanted as a function of frequency. 

2. _ Obtain a family of curves in values of frequency (one for 
each fixed frequency) for the magnitude of the output as function of the 
nonlinearity element values. This nonlinearity should be allowed to 
‘vary along a range that includes the values that are going to be read 
when it is expressed as a function of the magnitude of the voltage across 
its terminals. 


These curves may be obtained by hand computation or computer computa- 


+e 


i tion of the mangitude of the output function, for a certain input, using 


ef 


successive fixed values of the nonlinearity, which is treated as a para- 
meter. Using a computer program it is very practical to obtain these 
families of curves and the graphical method leads very simply to the fre- 
v quency response plot. 

3. Having obtained the "parameter versus magnitude" curves, 
it can be seen that for each fixed value of the parameter (which is rep- 
resented in.the graph by a line perpendicular to the parameter axis, and 
intersecting it at the chosen value), the magnitude values of the output 
function are given by reading the points along the magnitude axis corres- 


ponding to the intersection of the line representing the parameter value 


2e 


with each one of the curves of the frequency family. 
= 4, It was this concept that lead to the idea that the graph 


could be used: to read values of magnitude corresponding to the inter- 


a 





section of the frequency curves with any value of the parameter, in such 


a 19 








a way that if the parameter could be expressed as a function of the 
magnitude in a graphical manner (this for the case of nonlinear varia- 
tion of the parameter), and superimposed on the parameter plane curves, 
the frequency response of the network could be obtained by reading it 
directly from this graphical display. 

5. This is in fact what was done, and the present work in- 
cludes several types of networks with different nonlinearities to which 
the method was applied. At the same time the frequency response of these 
networks was obtained by means of simulation using a digital computer 
and the system IMB DSL 360. 

The results obtained, as was stated before, may be considered to 
fall within acceptable accuracy. 

6. The examples that follow show the method used in some net- 
works of general use, and the relationships between the parameter and 
the voltage across its terminals does not correspond, in some cases, to 
that most likely to be encountered (example, negative temperature co- 
efficient resistors). This was done to cover a large variety of func- 
tions, in the general sense, and as a way of illustration. 

Also, the cases considered have normalized values for components, 


as a means of simplification. 


3.2 Examples of Application of the Method. 
3.2.1 Second Order Network (Band-pass filter) 


The first case considered was the following network, (Fig. 3.1), 
in which the output resistor was a nonlinear function of the voltage 2 
across its terminals. 


The values used are normalized and the input voltage is considered 


to be a pure sinusoid of magnitude unity. 





Consider the network: 


R, L ¢ 
WI na 


é. P= Re ‘i a 


Pe Le 


Fig. 3.1 Band-pass Network 


Using normalized values for K , L and C it can be seen that its trans- 


fer function is: 


Te tol - Rz 
E-¢) (Rae Rj+sl+ Va 


or 


715) =z. € an ee 4 _ <5 


((, +Rics4Sloe- 1 Lest (Qt Ryeswa 


> 


ce: 
oe S24 (Arrar J 


where 
Rist, Leis C mj z 

One of these elements may be used as the constant parameter »« for 
PARAM 5 -_program purposes. 

Using equations (2-5) and (2-6): 

2 

w 8 ws 
@-ul Pe @+) > w™ 
This equation is the one that could be used to obtain the ‘magnitude 


mts Mtyjw) Tes w) = 


ce 4 


versus (2" curves, for successive values of frequency wW . In this work 


the PARAM 5 Program was used and the graph obtained shows (2 varying 


Td 


from 0 to 4.0 while the magnitude of the output voltage goes from 0 to 


8, and is shown on Fig. 3.2. 
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After the graph wan obtained, three cases of possible nonlinear 
variations of 6 with voltage were considered, and their graphical rep- 
resentation superimposed on the main graph. The values obtained by an 
ing the magnitude corresponding to each intersection of the fp lines 
with the frequency curves were plotted to obtain the magnitude-versus- 
frequency curve which in fact corresponds to the frequency response of 
the network. 

The same cases were simulated using digital simulation, and the 


results obtained for the frequency response are also shown. 


First Nonlinearity 


The first nonlinearity considered was of the type shown, 


(6) 


ae 


VeLT- MAGA TUDE 





Fig. 3.3 Beta as Nonlinear Function of Magnitude 
in which p is considered to maintain a constant value of 2, until the 
instantaneous magnitude of the voltage is 0.2, and after that to decrease 
constantly with a slope of -2/.3. 


The results obtained are tabulated and presented on the graph on 


Fig. 34. 





OUTPUT 
VOLTAGE 
MAGNITUDE 
-8 CURVE 1 Linear response 
CURVE 2 Nonlinear response 
\ Values checked by simulation 
6 
x A 
J V4 . 
| P r 6 CURVE § 
2 4 = 7 a 
CURVE 2 
Us 2, 4, 6. 6. 10. 


W RAD/ SEC 


Fig. 3.4 Frequency response of nonlinear band-pass network (CASE 1) 








+ 


Obtained "values for. frequency response of band pass network 


wt Output Magnitude ~ Output 
ee . PARAM 5 Magnitude 
(freq. Rad/Sec) Graphical Method Simulation 
1.0. _ “40 417 
2.0 370° .398 
3.0 - 335 373 


5.0 .280 | _ 300 





Second nonlinearity 


The second case considered, using the same network, was of the 


same type but with a different variation for.3°. 





(#) 
BETA 
Dr 
2. 
| 
lo . . 
: Vor. MAGMITUDE 
Fig. 3.5 Beta as nonlinear function of magnitude 
The slope now is a and the initial value for @® is 3.0. The results 


obtained follow, and the corresponding graph is included as Fig. 3.7. 








nate oe 


(freq. Rad/Sec) Output Output 











@ Voltage Magnitude Voltage Magnitude 
(PARAM 5): (Simu lat.) 
1.0 | 625 .618 
"90 7 875. - 565 
4.0: .450 "440 


5.0 405 | .400 
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Third Nonlinearity 


r 


of the filter with the nonlinear element increasing its value with higher 
values of voltage 
(2) 
BETA 
% 


2 


7?) re -6 VoLT. MAGMITUDE 


Fig. 3.6. Beta as Nonlinear Function of Magnitude 
The results obtained by simulation and graphical methods are as follows, 


and the corresponding graph is included as Fig. 3.8. 


Frequency Response of Band Pass Net, for Third Nonlinearity 





a Output a Output 
Freq. (Rad/Sec) Voltage Magnitude Voltage Magnitude 
(Simulation) (PARAM 5) 
1,2 ES 75 
2.0 713 70 
4.0 255 52 





The third case considered, with the same network, shows the behavior: | 
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ig CURVE 1 Linear response ; 
CURVE 2 Nonlinear response 
‘a Values checked by simulation 
6 
D oe 
_-» CURVE 4 
nN a: - 
~S 
CuRveE 2 
2 
0:0 - 2.0 4.0 6.0 8.0 10.0 
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W RAD/SEC 


Fig. 3.7 Frequency response of nonlinear band-pass network (CASE 2) 
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a Values checked by simulation 
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Fig. 3.8 Frequency response of nonlinear band-pass network (CASE 3) 


4, Use of the Method in the Shaping of. Network Responses. ma . 


One area that may be considered, of interest is the possibility of 


use of the present method and parameter plane curves to find what kind 


of nonlinear element should be introduced in a network in order to obtain. 


a certain frequency: response, different from that obtained with linear ~ 


elements in the network. 

a For a case like this, after obtaining the linear begquaniby response Z 
it may be seen from the respective graphical display chat a somewhat 
different curve shape is desired, corresponding bs different frequency 
response. This desired response may be cemndinprost on the parameter - 
plane, by plotting the chat buderyeaeuediiied. cacy ara luns that correspond 


_ Once in the parameter plane, these points define a curve for the 


parameter in terms of magnitude of the response,‘and may be.used to. ap- 
ei : ‘ gv 


of 


proximate a nonlinear function which would express the parameter in terms 


. . . , , 4 . ‘ 
of magnitude. Such a nonlinear element,if inserted in the network, would 


provide the frequency response which originated the same functional re- 


lationship. 


The degree of accuracy and the closeness between the desired values 
- oad 
and those obtained depend on the approximations made when expressing the 
parameter in terms of magnitude in order ‘to make it correspond to the 


series of points transported from the frequency response plane. 


This method will be illustrated with the following examples: 


a 


4,1 Reshaping of frequency response curve of a. band-pass network, by 


i introduction “of -a nonlinearity. ve 
Consider. the network used {n example: 3.3.1; the frequency response * 


for this network is easily obtained by the graphical method explained 








previously and is shown in the graph in Fig. 4.2, for linear values of 
the elements. In this case it was considered that the output resistor 
had the invariant value of 2.0 units, and the curve obtained is refer- 


red to as number l. > ae 


Fig. 4.1 Network Considered in the Present Example 
If it is desired to sharpen the response of the present network to a 
different cvs as for example the one indicated in the graph of Fig. 
4.2, curve 2, the points may be taken from the new curve and transported 
to the parameter plane, and plotted there. 

This was done in the graph on Fig. 4.2, as indicated by curve n@ 2, 
and the points corresponding to the desired response were plotted in 
curve n+. 4 on the graph on Fig. 3.2, which refers to the network being 
considered. These points are indicated by the symbol E] ° 

‘The relationship between magnitude of voltage and values of non- 
linearity were read from this same graph, 

This having been done, it was decided to use the linear, time-wariant 
resistor Ry as a time-variant element, defined by the following equation 


of a straight line in the parameter plane: 


Ry = (output voltage magnitude) x 2.22 + .5 (4-1) 


ro 


It may be seen that if the output resistor really behaves in accordance 


with equation 4-1, the frequency response is the desired curve, when:-— 


transported back to the frequency-response graph. 
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CURVE 1 Linear response (8 = 2, 
CURVE 2 Desired nonlinear response 


A Values checked by simulation 
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-The present example was checked by the use of digital simulation, 
and by verifying that the frequency response obtained by simulation of 
a nonlinear network, with the nonlinear element expressed by equation 4-1, 
was a reasonable approximation to the desired frequency response. | 


For this’ case, the following results were obtained: 





Linear Frequency Response Ry = 2.0 Units 
“w = 1. tnd/eed Magn = .19 W = 4. rad/sec Magn = .42 
w= 5 Magn = .6 w= 5.0 Magn = .35 
w= 1.0 Magn = .67 (v= 6.0 Magn = .305 
We 2,0 _ Magn = .6 _ We 7.0 Magn = .27 
5.55 ws 8.0 Magn = +24 


‘W= 3.0 Magn 








After reshaping the response, the following values were obtained 


Desired Frequency Response for Nonlinear Network. 





Magnit. Obtained 


: »b Simulation : 


WRad/Sec Desired Magnitude Corresp. R 








; 4.2 Frequency Response and its Reshaping, for a low-pass Filter with Load 


As another illustration consider the following network. © - 


R be R , L OSes oe 
t A a 
4 f , * 


Fig. 4.3 Second Order Low-Pass Filter. 





with two possible variable parameters, ™ and ke , or with a variable 


parameter 6 which could be used in conjunetion with a different load 


& for each case. 


The transfer function is obtained as: 





Ns) = Eo = wf 


~~ t CR*K B+ 2LER+ Leo] +s [oucrup’ritc] +s# falc] 


é 
i 


(4-2) _ 


It is desired to analyze the transfer function behavior, or the 


frequency response for a particular value of o& when ( is a nonlinear 


function of the output voltage. This can be done by the PARAM 5 method 


previously explained, and the following. results were obtained for the 


simplified case, with all values normalized, in which 


a? 


a & ai BR el 


‘The transfer function becomes: 


na oe « ae + 
at «+s[ 20g + ara ]+ SBupr2+« | +s [z«p+i | +s# [« é | : _ = 


sa We 
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and for & = 1, this reduces to: 


Is = 3+ sfea+4 | +5” [3e+3] esi [apes +s*[ py 


(4-4) 
Using Eq. 4-4 for the parameter plane program the curves relating and 
magnitude of output voltage were obtained. Two cases were considered, 
one for each type of nonlinear behavior of eB » as shown on the graph 


on Fig. 4.4. 


. 4.2.1 For the first case consider & as a function of magnitude of 


output.voltage, as given by the following relationship: 


(8) 
Be7A 4. 





VOLT. MAGNITUDE 
Fig. 4.5 Beta As Nonlinear Function of Voltage 
mi Magnitude 
Beta (8) =4 (or output voltage magnitude <= my) 


or Beta ( f&) = -9.51 x (output voltage magnitude) + 4.951 
Employing the usual procedure of superimposing this graphical relation- 
ship on the garemater plane curves, touree 1) the frequency response was 
obtained and the results checked by simulation of the etapa, with the. 


graphical presentation included as Fig. 4.6, ‘ 
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VOLTAGE 
MAGNITUDE 
3 CASE 1 as nonlinear function of voltage 
f magnitude 
ry Values checked by simulation 
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Fig. 4.6 Frequency response of a two-section low-pass 
Network 
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ov 


on Fig. 4.4, (Curve 2). A 


Results. Obtained for Freq. Response of.a 2-section Low-pass Network 





Case 1 
W Rad/sec. Out. Volt. Magnitude -. Out. Volt. Magnitude « - 
(Parameter Plane) i (Simulation) — ee 

2 325 sy . | - 
a. | .158 0 a 

1.0 120 119 

Le * ,095 -0961 — 

2.0 .02 .0207 


« 


The above results are plotted in the graph on Fig. 4,6,-and are considered 


“to be a reasonably good approximation. 


“4.2.2 


-— wae a second illustration consider a different type of nonlinearity 


‘for B - . ae . . : ‘ 





AEF ae il 


co 4 “495 Vor. MAGHITYDE 


Fig. 4.7 Beta-as Nonlinear Function of Magnitude 
As shown above, 8 is defined by ‘the equations: 


{s = .17 for output voltage magnitude “¢-4 
" 5 = .i¢ x (out. voltage magn.) for .1< output voltage. 
Magnitude <£. .25 


Bp = .4 for out. voltage magnitude om -.25 


Using the usual method the following results. were obtained from the. graph 


2° * 





Results obtained for Freq. Response of a two-section, low-pass network 


Case 2 
wu) rad/sec Out.Volt. Magnitude Out.Volt. Magnitude 
Peart oe (Pare Blane) eG imutationy 
2 2295 313 
8 - 180 -208 
1.0 .170 174 
1.2 . 150 . 160 
2.0 .050 .050 
2.4 | .020 -023 


These results are displayed on the following graph, on Fig. 4.8. 


4.2.3 Reshaping of Response of Networks. 

Now consider the possibility of using the method as a means of 
changing the response of the network by dnewoduct ten of a nonlinearity, 
as was explained in Section 4.1, 

In this example the linear response of the network was obtained for 
fp = 2.4, from the parameter plane curves, and it was considered desir- 
able to flatten the response in a way that would give constant magnitude 
’ response up to a frequency of “= .8 rad/sec, as indicated by Curve No, 2 

in the graph on Fig. 4.9. The points read from this curve were eee 
ported*to the parameter plane graph and plotted as Ch) ote was chosen to 
approximate them by two segments of straight line as shown by curve 3. 
The resulting frequency response for the approximation of P as a function 
of output voltage was plotted in the graph »n Fig. 4.9, indicating that 
the network had in pak a response reasonably close to what was desired. 


The defined relationship between p and output voltage was simulated 


and the results obtained are shown in the following table. 
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Fig. 4.8 Frequency response of a two-section low-pass network 








CURVE 1 Linear response © = 2.4 
ovTPpuT 


VOLTAGE CURVE 2 Nonlinear response from parameter-plane 
MAGW ITU DE CURVE 3 Desired shape for response curve 
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A Values checked by simulation 
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Fig. 4.9 Frequency response of two-section low-pass network 
with load 
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4 ee Results Obtained in Reshap ing the Frequency Response of: a Section Low-Pass Filter ! 
W rad/sec Desired Magnitude ~ “" " Aet.’ Introd... ~~ Expected Response Obtained Response 






Response Nalue of from Param. . Plane 
2 : 5 . fags 335 
wo 
AS | 
_ £325 5:5 330 350 
3 oer 
= & 
.8 = x a 30 330 
1.0 £275 + & 2275 .280 
. s : : 
alist ee Le * 2 «252. «252 
Ee as : r j ; = 
1.4 : .225 ve @ - 225 i 
1.8 .150 aD .150 150 
2.0 .120 3 .120 ak 





Di. Frequency Response of Networks with Output Voltage not Directly 


Related to the Nonlinear Element. 

As an extension of the present work, the case can be considered where 
the nonlinear element is not directly across the output acinade abut may 
instead be in any other place in the network. 

It can be seen that in a case of this kind, the method cannot be ap- 
plied directly, but instead can be used to determine the values that the 
nonlinear element will take on as a function of the input voltage. Then 
using these values for a point-by-point solution of the main transfer 
function of ultimate interest, the frequency response of the network can 
be determined. 

To do this, a secondary transfer function is defined, relating the 
input voltage and the voltage across the nonlinearity, and is used in a 
parameter plane context to obtain the nonlinearity as a function of 
frequency (and consequently of the voltage magnitude across itself). 

This is accomplished by superposition of the graphical relationship be- 
tween the nonlinear element and the voltage across it on the parameter-~- 
plane frequency curves, as was done before. 

With this relation between frequency and the value of the nonlinear 
element, it is possible to enter the main transfer function relating in- 
put voltage to output voltage, and compute output voltage for each fre- 


quency. 


5.1 Indirect Nonlinear Frequency Response. 
As an illustration, consider the network for which C is considered 


R 








Fig. 5.1 Network for Indirect Response 
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as the nonlinear element, i) , 


The flow graph for this network is 
| Be 


ESF Lad 
1% Ae 


4 
vA 
ee eee 


Fig. 5.3 Flow Graph for Network of Fig. 5.1 


A 


and it can be seen that 


Tis) - Mats . Pee 5-1 
6 Es) (ReRjse + 4 on 


or 


7 (9 a , t (5-2) 
(Ri + Ri)s@ + 4 


The parameter-plane curves for this transfer function T, (s) were 
obtained (Fig. 5.2), and the nonlinear relationship between f and 
magnitude of voltage across it, for this particular example, was super- 
imposed on these curves, making it possible to obtain successive values 


of (es as a function of frequency. 


From the flow graph for the network, it is seen that 


(5-3) 


[) - Ent) _ SCR, - S$ R, 


Evi) (Rat@)sc +s (Rr Ropar! 





a@ 


This is the principal transfer function, giving output voltage in terms 
of input voltage. 


It is also seen that 


E, (5-4) 
ie tk)” a al | 
Equation 5-4 represents the frequency response of the network, giving 
magnitude of output voltage in terms of frequency. For the present case, 
(3 is nonlinear, but also has been obtained as a function of frequency, 
from the parameter plane. . 
Using this result, it is possible to obtain the frequency résponse 
of the. network, by hand calculation for several values of W . 
In this example, f was chosen as the following function of Ve 


(A) 
BETA 





MAGNITUDE OF Vp = 


Fig. 5.4 Beta as Nonlinear Function of ‘Magnitude 
B= -2 for magnitude of Ve = .2 
(= -l + .5 x (magnitude of Vp ) for magnitude of Wr .2 
and R, * IF R, = 3 and JE: | = 1. 
The results obtained are summarized in the following table, and it 


was considered that the simulation results agree, within acceptable range, 


with the values from the method presented. 
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Results Obtained by Using Parameter Plane Method on Indirect Nonlinearity, 
for Frequency Response of First Order Filter 





w rad/sec “VeMagn. Ve Magn. B From l&] From | |Eo! From 
From Param Plane (Simulation) From Param 5 Formula (5-4 Simulation 
1. -575 -65 38 .63 +55 
69 -655 
.725 -716 
.73 
74 74 
-745 744 


-748 747 





a®@ 


(Oe Conclusions. 
In the present work it was shown that it is possible to obtain the 
frequency response of networks with one nonlinear element, (whether or 


not ‘it is directly across the output voltage) by ‘a graphical method using 


the parameter plane concept: Also, the fact that computer programs were 


avai lable, Mae 3 easy to obtain the parameter plane. ee thereby 
avoidtg the. ‘somewhat mote extensive calculations inservbay describing- 
function techniques. ‘ Pe 

The method may be considered reasonably’ accurate, giving acceptable 
results accurate within the usual approximations. It was also shown that 
it provides a rapid way of defining the kind of nonlinear element needed 
to reshape the network response in a desired way, the only limitation here” 
being the practical realization of the nonlinear element found to be re- 
quired. 

As areas for extension of this work it would be of interest to study 
the applications of the parameter-plane curves in connection with the 
jump resonance phenomena, for nonlinear networks, and the possibility of 
modifying the response by changing the nonlinearity so as to avoid three- 
point intersections with the frequency curves in the parameter plane, 
Also, the investigation of the behavior of networks and the applicability 
of the method for cases where the input voltage has magnitude different 
than unity, that is, large or small values of voltage. The usefulness 
of the method for cases where there is more than one nonlinear element 


present, and the combination of their effects for obtaining a certain 


response, is another subject which may prove to be of interest. 
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